Pump-induced Exceptional Points in Lasers 
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We demonstrate that the above-threshold behavior of a laser can be strongly affected by excep- 
tional points which are induced by pumping the laser nonuniformly. At these singularities, the 
eigenstates of the non-Hermitian operator which describes the lasing modes coalesce. In their vicin- 
ity, the laser may turn off even when the overall pump power deposited in the system is increased. 
Such signatures of a pump- induced exceptional point can be experimentally probed with coupled 
ridge or microdisk lasers. 

PACS numbers: 42.55.Sa, 42. 55. Ah, 42.25.Bs 



The interest in physical systems described by non- 
Hermitian operators has recently been revived by a num- 
ber of pioneering experiments [1-7] which demonstrated 
the rich physics induced by the system's non-Hermiticity. 
In this context exceptional points (EPs), at which non- 
Hermitian operators are defective and pairs of their eigen- 
states coalesce [8], are of particular interest. In the vicin- 
ity of these EPs the eigenvalues display an intricate topol- 
ogy of intersecting Riemann sheets, leading to nontrivial 
geometric phases when encircling the EP parametrically, 
as well as to problems with mode labeling [4-11]. 

Photonic structures in the presence of gain or loss are a 
natural arena in which EPs can play a role, since they are 
described by a non-Hermitian Maxwell operator, arising 
from a complex dielectric function. Striking examples 
which have received much attention recently are optical 
systems with balanced gain and loss, for which this oper- 
ator is invariant under parity (V) and time-reversal (T) 
[12]. In such systems, as the real wave vector, k, of the 
input radiation is varied, the electromagnetic scattering 
matrix undergoes a transition at an EP from a "flux- 
conserving" behavior to an amplifying/attenuating be- 
havior [13]. This transition leads to a number of novel op- 
tical properties [1, 2], such as power oscillations of light, 
beam steering, and self-sustained radiation [14-16]. 

Lasers are prototype non-Hermitian systems, with the 
additional complexity of being described by a nonlinear 
wave equation above their first oscillation threshold. Un- 
til very recently, however, no connection had been made 
between pumping a laser and the phenomena associated 
with EPs. The reason for this, as we shall see below, is 
that the non-Hermitian operators describing lasers do not 
reveal EPs when the laser is pumped uniformly in space. 
Two previous works have hinted at EP physics in non- 
uniformly pumped lasers. First, a work on VT optical 
systems has shown that a cavity with balanced regions 
of strong gain and strong loss can actually lase after the 
PT-breaking EP has been crossed [13, 17]. However, 
the properties of this unusual laser (which has yet to be 



demonstrated experimentally) are not influenced by the 
EP above threshold. Another laser system which exhibits 
EPs is composed of very- low Q cavities, which can form 
novel surface-state lasing modes [18]. Again, however, 
these modes are created well away from the EP and their 
lasing properties are not affected by it, nor is this a prac- 
tical system for experiments. 

In the current Letter, we propose a class of non- 
uniformly pumped laser systems, which are achievable in 
practical laser structures and which show the dramatic 
role of an EP above the laser threshold in controlling and 
even switching off the coherent emission. The systems 
we study consist of two coherently coupled laser cavities 
as realized, e.g., by two coupled ridge lasers (see insets 
in Fig. 1). The cavities are assumed to have large loss in 
the absence of pumping. Their unusual behavior when 
pumped is shown in Fig. 1 (red curve) [19]. First the 
left cavity is pumped (0<d< 1) until it begins lasing in 
an asymmetric mode which emits predominantly to the 
left; then, pumping is added to the right cavity and in- 
creased (1 <d< 2), which, surprisingly, causes the laser 
first to shut off and then to turn on again as a symmet- 
ric mode (details of the pumping scheme are given in 
the caption). This behavior occurs for a specific value of 
the center and width of the gain curve. Quite different 
behavior occurs if the gain center is shifted by a small 
amount compared to the free spectral range. If shifted 
down (yellow curve), the system never lases at all in the 
asymmetric mode; if shifted up (blue curve), it does not 
ever switch off, but the mode slowly evolves from asym- 
metric to symmetric with a strong intensity dip after the 
initial turn-on. We will show that this counter-intuitive 
dependence on the pump strength and high sensitivity to 
the gain curve is completely controlled by an EP in the 
relevant non-Hermitian operator for the laser problem. 

The above results have been obtained from the sta- 
tionary solutions of the nonlinear Maxwell- B loch equa- 
tions, which are the key equations in semiclassical laser 
theory [20, 21]. To calculate these solutions, we em- 
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FIG. 1. (color online) Intensity output of a laser system con- 
sisting of two ID coupled ridge lasers, each of length 100 /am 
with an air gap of size 10 /am and an (unpumped) index of re- 
fraction n — 3 + 0.13i. For < d< 1, the pump in the left ridge 
is linearly increased in the range < D < 1.2, and for 1 < d < 2, 
the same is done in the right ridge while keeping the left ridge 
at D = 1.2. The points show the results determined from a 
finite- difference time-domain solution of the Maxwell-Bloch 
equations (using 7y = 0.02 mm -1 ). The solid lines are the 
solutions as given by SALT (see the text). The three differ- 
ent colors (red, blue and yellow [19]) represent the results for 
different gain curves with a width of 2j± = 4 mm -1 centered 
around k a = {92.4, 94.6, 96.3} mm -1 , respectively. Although 
the overall pump power increases monotonically over the en- 
tire range of the laser is found to turn off in the range 
1.55 < d < 1.68 due to the influence of an exceptional point. 
The insets show the electric field amplitude in the laser at 
the pump steps d — {1,2} for the gain curve centered around 
k a = 94.6 mm- 1 . 



ployed the newly developed steady-state ab initio laser 
theory (SALT) [22-24], which is exact in the single- mode 
regime of interest here (except for the standard rotat- 
ing wave approximation). Because of the novelty of the 
behavior studied here, we independently confirmed the 
SALT results (see lines in Fig. 1) using the more famil- 
iar finite-difference time-domain method (see points in 
Fig. 1) [25]. The observed reentrant behavior, however, 
cannot be understood in the time domain, requiring, in- 
stead, a frequency-domain approach like SALT. In ID 
and 2D, the SALT equations are coupled scalar nonlin- 
ear wave equations for the lasing modes ^(x) and their 
corresponding lasing frequencies of the form [22] 



V 2 + ^£(x,{^,^}) 



^(x)=0. 



(1) 



Here the complex dielectric function s(x, {k Vl = 
£ c (x) +^(x, {k Vl ^ v }) comprises a contribution from the 
passive cavity e c as well as a gain contribution e g induced 
by the external pump. This active part 
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contains modal interactions and gain saturation, due to 
spatial hole burning through the nonlinearity given by 
[1 + Ysv T A^A 2 \~ l i where T u = T(k v ) is a Lorentzian 



function of width 27^, centered at the atomic transi- 
tion frequency k a and evaluated at the lasing frequency 
k^ (c = 1). Crucially, e g in Eq. (2) also contains the 
spatially varying inversion density £>o(x, <i), the "pump" 
[26], where we have introduced an adiabatic parame- 
ter d. This leads us to the concept of a "pump tra- 
jectory" (see, e.g., Fig. 1) consisting of a sequence of 
adiabatic variations in both the pump profile and ampli- 
tude, parametrized by d and satisfying <C 7y to en- 
sure adiabaticity (7N is the longitudinal relaxation rate) 
[20, 25]. We obtain the solutions to Eq. (1) along such a 
pump trajectory for both the laser frequency k^ and the 
corresponding electric field as the fixed points of a 
self-consistent iteration carried out for each given value 
of the pump Do( x , <i) [27]. With this procedure, we find 
the full nonlinear solutions for the emission intensities 
above threshold; as we now demonstrate, the qualitative 
behavior shown in Fig. 1 is, however, determined already 
by the linear threshold conditions for a given pump tra- 
jectory. 

To show this, we introduce a generalization of the 
threshold equations of SALT, which are posed in two 
steps. First we define a non-Hermitian eigenvalue prob- 
lem, 



V 2 +& 2 £ c (&)+7M(M)A)(x,d) 



>(M;x) = 0, (3) 



where r] n (k,d) is the complex eigenvalue [22, 27], and 
r] n (k,d)Do(-x,d) plays the role of e g in Eq. (2) at thresh- 
old, where the nonlinearity vanishes. Equation (3), which 
is solved using non-Hermitian "constant flux" boundary 
conditions at the laser boundary (d x u n = ±iku n ), is a 
generalization of the threshold constant flux (TCF) equa- 
tion [22] for pump trajectories; the {u n ,r] n } are the TCF 
states and eigenvalues. 

The first laser mode is determined out of the countably 
infinite set of r] n (k,d) as the one which satisfies 
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k a + iji. 



(4) 



for the smallest value of <i M , which, from Eq. (4), also 
yields the threshold lasing frequency k^. The lasing 
mode ^(x) at the threshold is then just given by the 
corresponding eigenstate u n (k^, d^\ x). The linear non- 
Hermitian TCF equation, Eq. (3), can, however, be 
solved for arbitrary values of k and d. We are thus natu- 
rally provided with the two continuous parameters neces- 
sary to demonstrate the presence of EPs in the complex 
spectrum of the TCF operator [8]. This will bring us 
to the central insight that the TCF operator which de- 
termines the lasing thresholds contains the EPs which 
control the lasing behavior shown in Fig. 1. 

To show this explicitly, we rewrite the complex Eq. (4) 
as two real conditions which neatly separate the role of 
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the TCF spectrum and the gain curve 
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lm(r] n ) 0, 



Im(ry n ) 



7± 



Re(r/ n ) 



(5) 



The left equation has not been analyzed in previous work 
on SALT; it has the remarkable property of being en- 
tirely independent of the gain curve parameters k a and 
7^ (which do not appear in the TCF equation). Hence, 
it defines a threshold boundary in the /c, d plane divid- 
ing lasing and nonlasing regions, solely from the values 
of {r] n }. To conveniently capture the entire k,d land- 
scape of possible lasing solutions, we define a function 
/(fc, d) = min n [ \r] n \ 2 + Im(r/ n ) ], which has the property 
that the contour f(k,d) = is the locus of all possible 
lasing thresholds. The actual lasing thresholds on this 
contour are then determined by the gain curve param- 
eters in the right-hand condition of Eq. (5). In Fig. 2, 
we show the relevant regions of f(k,d) above (/ < 0, 
green, within black contour) and below (/ > 0, red, out- 
side black contour) threshold, as well as the threshold 
boundary (/ = 0, black contour). When / < 0, the TCF 
equation does not give the lasing solutions and full non- 
linear SALT is needed to find the intensity and frequency, 
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FIG. 2. (color online) Contour plot of the function f(k,d), 
which indicates the parameter regions of k and d where the 
laser in Fig. 1 is above threshold (green, inside black con- 
tour) or below threshold (red, outside black contour). The 
threshold condition, /(&, d) = 0, is satisfied at the solid black 
contour. The EP in the center of the plot pulls the laser below 
the threshold in its vicinity. The frequency dependence of the 
solutions of the nonlinear SALT equations (dashed lines) and 
the corresponding laser thresholds (dots on the solid black 
contour) are provided for the same gain curves (shown in the 
top panel) as in Fig. 1. The evaluation of /(/c, d) involves two 
interconnected Riemann sheets (see Fig. 3), resulting in a cut 
between the two sheets (see the white line right above the 
EP). The two black dots at the upper edge of the plot repre- 
sent the resonance frequencies of the passive cavity system. 



as was done for the curves in Figs. 1 and 2. 

The most striking feature of the contour plot shown 
in Fig. 2 is the existence of a substantial inclusion of 
subthreshold region (red, outside black contour) in the 
midst of the superthreshold region (green, within black 
contour). At the center of this inclusion is a local maxi- 
mum of /; this point corresponds to an EP of the TCF 
operator. For the middle gain curve (red), there are three 
threshold solutions, two very near the EP. If we order 
these solutions according to increasing pump, <i, we find 
that df/dd < at the first and third thresholds and 
df/dd > at the second. As a function of d, the laser 
thus turns on-off-on, just as found in Fig. 1. Hence, non- 
linear effects play no role in the qualitative behavior of 
this laser; they only determine the amplitude of the lasing 
emission in the green region. All the interesting behavior 
is controlled by the EP: It completely suppresses lasing 
in its vicinity of k and d-parameters, thereby causing the 
reentrant lasing (red curve) as well as the nonmonotonic 
but continuous lasing emission (blue curve), which is in- 
fluenced by its proximity to the EP. 

In the entire k and d range shown in Fig. 2, exactly only 
two eigenvalues rji^ contribute to f(k,d). In the coupled 
cavities studied here, these two eigenvalues are associ- 
ated with a symmetric-antisymmetric doublet of passive 
cavity resonances (see black dots in Fig. 2). The dou- 
blets are separated from each other by approximately 
the free spectral range of a single cavity. If both cav- 
ities are uniformly pumped with the same gain curves 
as in Fig. 2, the laser emits in modes each of which can 
be associated with one of the passive resonances (not 
shown). However, no EP and hence no pump-induced 
inclusion or nonmonotonic laser emission appears. Only 
when nonuniform pumping is applied may a pair of eigen- 
values coalesce at an EP. 

In Fig. 3, we plot the k and d dependence of the real 
and imaginary part of the two eigenvalues 771,2- Focus- 
ing on the parameter region around d— 1.55, where the 
laser turns off, we find that the parametric dependence of 
the eigenvalues on k and d shows the typical topological 
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FIG. 3. (color online) (a) Real and (b) imaginary parts of the 
two eigenvalues 771,2(^,0?) which are closest to the threshold 
in the k and d parameter region around the laser turn-off 
observed in Fig. 1. The eigenvalue surfaces display the typical 
structure of intersecting Riemann surfaces, centered around 
an EP. 
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FIG. 4. (color online) Laser light intensity emitted from two 
coupled microdisks (see insets) as a function of pump param- 
eter d. With the same pumping scheme as in Fig. 1 (maxi- 
mum pump in each disk: D = 0.71) the laser is also found 
here to turn off for increasing d. The disks have a radius 
of 45 /im, an inter cavity spacing of 10 \ivn and an index of 
refraction n = 3.67 + 0.09i. The gain curve parameters are 
k a ~ 2.82 THz and j± ~ 0.29 THz. The insets show the ab- 
solute value of the electric field of the green (dashed) lasing 
mode at the pump points d = 1 and d = 2. 



structure of two intersecting Riemann sheets around an 
EP located at k « 95 mm -1 , d ~ 1.6. Away from the 
EP, the laser can always "choose" the lower threshold 
solution for lasing (upper part of the Riemann sheet), 
but at the EP only one solution is available, which is a 
compromise between the low and high threshold solution. 
This gives rise to a local minimum in the effective gain 
and inhibits lasing in the vicinity of the EP. In this way, 
the essential physics of an EP directly translates into the 
physical effects observed in Fig. 1. 

To confirm that such a behavior is generic for coupled 
microlasers, we consider a device consisting of two cou- 
pled circular microdisks [28, 29] which has recently been 
implemented as an electrically controllable "photonic- 
molecule laser" [30]. For this more challenging compu- 
tation, we have incorporated in the SALT an advanced 
higher-order finite element method discretization tech- 
nique [31]. When applying the same pumping scheme 
as in Fig. 1 to such a coupled microdisk laser (see in- 
sets of Fig. 4), we again find that the laser may turn 
off when increasing the overall pump power (see main 
panel of Fig. 4). Note, however, that more than one 
mode is lasing both below and above the turn-off due to 
a quasidegeneracy of modes. This example demonstrates 
that the effect of EPs on the above-threshold character- 
istics of a laser should be generally observable in coupled 
microlasers. 

The success of a corresponding experiment hinges on 
striking the right balance between the coupling strength 
and the amount of gain or loss in each disk. Gain or loss 
is fixed by the cavity material as well as by the applied 
pumping strength. Microlasers built from active quan- 
tum cascade structures (QCLs) seem to be ideally suited 
to test this effect, since they can be well-controlled in 



their shape and hence in their coupling. With the active 
medium being situated within a lossy metal waveguide 
[32], a central requirement for observing this effect, i.e., 
that the absorption length in the passive cavity is smaller 
than the round trip length, has already been experimen- 
tally realized for a variant of such lasers [33]. Also, a 
sufficiently strong inter-cavity coupling has been demon- 
strated for such coupled QCLs [30, 33], suggesting that 
an experimental realization might be within reach. Note 
that we have used realistic material parameters from the 
experiment [30] in the employed MB equations which typ- 
ically describe QCLs very well [34]. 

To summarize, we have shown that the presence of 
an exceptional point in the lasing equations for coupled 
microlasers may be "brought to light" by a suitable vari- 
ation of the applied pumping. We find that, in the vicin- 
ity of the exceptional point, the output intensity of the 
emitted laser light is reduced, although the applied pump 
power is increased. This unorthodox and robust lasing 
effect may even turn off a laser above its lowest thresh- 
old. We suggest an experimental realization of this effect, 
which would open up many more possibilities to study 
the rich physics associated with EPs in the pump depen- 
dence of a laser. 
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